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Abstract. In earlier work, Jockusch, Propp, and Shor proved a theorem de- 
scribing the Umiting shape of the boundary between the uniformly tiled corners 
of a random tiling of an Aztec diamond and the more unpredictable 'temperate 
zone' in the interior of the region. The so-called arctic circle theorem made 
precise a phenomenon observed in random tilings of large Aztec diamonds. 

Here we examine a related combinatorial model called groves. Created by 
Carroll and Speyer as combinatorial interpretations for Laurent polynomials 
given by the cube recurrence, groves have observable frozen regions which 
we describe precisely via asymptotic analysis of a generating function. Our 
approach also provides another way to prove the arctic circle theorem for 
Aztec diamonds. 



1. Introduction 

Groves came into existence as combinatorial interpretations of the rational func- 
tions generated by the cube recurrence: 



fi.j,kfi — l,j~l.k~l — fi~lj.kfij — l.k~l fi.j — l,kfi — lj.k — l ~^ fij^k—lfi — lj- 



l.fc, 



where some initial fmictions are specified. Typically, /ij,fc ■— Xi^^k for some choice 
of {i,j,k) G Z'^ called the initial conditions. Fomin and Zelevinsky 0] were able 
to show that for initial conditions satisfying some basic requirements, the rational 
functions generated by the cube recurrence are in fact Laurent polynomials in the 
Xi^j^k- The introduction of groves by Carroll and Speyer P gave a combinatorial 
proof of the surprising fact that each term of these polynomials has coefficient 
The main results in this paper only apply to the family of groves on standard initial 
conditions as described in Section ll. 11 ^ 

Before getting into the details of groves, let us first describe the motivation for 
this paper: random domino tilings of large Aztec diamonds. An Aztec diamond of 
order n consists of the union of all unit squares with integer vertices contained in 
the planar region {{x, y) : \x\ + \y\ < n + 1}. A domino tiling of an Aztec diamond 
is an arrangement of 2 x 1 rectangles, or dominoes, that cover the diamond without 
any overlapping. A random domino tiling of a large Aztec diamond consists of 
two qualitatively different regions.^ As seen in the random tiling in Figure ^ the 
dominoes in the corners of the diamond are frozen in a brickwork pattern, whereas 
the dominoes in the interior have a more random, temperate behavior. It was shown 
in and that asymptotically, the boundary between the frozen and temperate 
regions in a random tiling is given by the circle inscribed in the Aztec diamond. 



^Herein we will invoke some of the basic properties of groves without proof. For such arguments, 
as well as a general treatment of groves and the cube recurrence, the reader is referred to 

^By random we mean selected from the uniform distribution on all tilings of an Aztec diamond 
of order n, though other probability distributions may be considered as well. See |2]- 
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Since everything outside the circle is expected to be frozen, it is referred to as the 
arctic circle. 




Figure 1. A random domino tiling of an Aztec diamond of order 64. 

In this paper we shall see that groves on standard initial conditions exhibit a 
very similar behavior. A grove, however, is not a type of tiling. As the name may 
suggest, a grove is in fact a collection of trees. From our point of view, groves 
are spanning forests on a finite triangular lattice satisfying certain connectivity 
conditions on the boundary. We will show that outside of the circle inscribed in 
the triangle, the trees of a large random grove line up uniformly. 

Despite their superficial differences, groves and random domino tilings of Aztec 
diamonds are linked by more than their asymptotic behavior. In fact it seems that 
their asymptotic behavior is similar because they share a deeper link. The paper of 
Carroll and Speyer U establishes that groves are encoded in the terms of a Laurent 
polynomial given by the cube recurrence. There is a more general form of the cube 
recurrence: 

fi,j,kfi-l,j-l,k-l — CK/'i-l ,j,fc/ij-l,fe-l + /3/ij-l,fc/i-l j",fc-l + 1 fi,j M-1 fi-l,j -l,k 

where a, (3, 7 are constants, li a — (3 — j — 1 we have the original form of the 
cube recurrence from whence come groves. If a — /? = 1 and 7 = 0, we have (after 
re-indexing), the octahedron recurrence: 
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with which we may encode tihngs of Aztec diamonds. In Section |21 we wiU show 
how the polynomial (?o,o,n yields all tilings of an Aztec diamond of order n and we 
will describe the role that this recurrence plays in the large scale behavior of such 
tilings. 

While the octahedron recurrence is important to us, it has not been extant in the 
study of tilings of Aztec diamonds in the past. Rather, a local move called domino 
shuffling has been used. Domino shuffling was introduced in ,3.^ and is generalized 
in It provides a method for generating tilings of successively larger Aztec 

diamonds uniformly at random, and has been at least implicit in all probabilistic 
analysis done to date. Section fl.3l will introduce an analogous local move for groves 
that we call grove shuffling. Like domino shuffling, it will be key to our analysis. 

For each of the two models discussed we have a global perspective and a local 
perspective. Laurent polynomials tell the global story: all groves are encapsulated 
in ^0,0.0 (from the cube recurrence), all tilings in (7o,o.n 

(from the octahedron re- 
currence). A specified shuffling algorithm tells the local story. In this paper we 
combine these two points of view to build generating functions (for tilings of Aztec 
diamonds as well as for groves), with which we can study asymptotic behavior. 




Figure 2. A portion of Q{b). 
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1.1. Groves on standard initial conditions. The standard initial conditions 
of order n specify a vertex set I{n) = C{n) U B{n) where C{n) — {{i,j,k) e Z"^ | 
-n-1 < i+j + k < -n+l,i,j,k< 0} andB{n) = {{i,j,k) e | i+j + k < -n-1; 
hji k < 0; and or k — 0}. We draw its projection onto the plane R'^/(l, 1, 1) 
as shown in Figure [5] for the case n = 5. One way to generate all groves of order 
n is to set fij^k ■= Xij^k for all {i,i,k) g 1{n), and compute /o,o,o- Each term 
in the resulting Laurent polynomial defines a grove as follows. Let Q{n) be the 
graph on the vertex set T{rL) where vertex (i, j, k) has as its neighbors the vertices 
I{n) n {(i ± l,j ± l,fc),(i± 1, j,/c± l),(i,j± l,fc± 1)}. Pictorially, edges oig{n) 
connect vertices that lie diagonally across a rhombus. In Figure |2 the graph g{b) 
is made up of the lighter edges and the dark vertices. 

As established in 1 , the terms in /o.o.o ai'e Laurent monomials of the form 



m{g) 



^dog(i,j\fc)-2 



where deg(i,j, fc) g {1,2,. 
We have the following 



(ij\fc)el(n) 

6} is the number of edges connected to vertex (i, j, k). 



Definition 1. The grove g defined by m{g) is the unique subgraph oj Q{n) con- 
taining no crossing edges such that vertex {i,j,k) in X{n) has exactly deg{i,j,k) 
incident edges. 

The uniqueness of the grove determined by each monomial is a consequence of 
Theorem 3 in For example, /o,o.o on 1(2) is 

-1,03^0,0,-1 ^ a^-i,o,-ia;o,-i,o ^ a:o^_i^_ix_i^o,o 
2^-1,-1,-1 x-i-i^-i a;_i,-i,-i 

and the corresponding groves are shown in Figure |21 






Figure 3. The three groves of order 2. 



For a more interesting example, one term of /o,o.o on 1(5) is 

a;-3,0, -23^-2, -l,-ia;-l, -3, 03^0,-2, -2 
a;-3, -1,-22^-2, -3, -12^-1,-2,-2 

Its corresponding grove, g, is shown in Figure^ We can observe some connectivity 
properties of this grove that in fact hold for all groves. Every vertex on the boundary 
of C(n) (where cubes have been pushed down) is connected to another vertex on 
the boundary of C(n) if and only if those vertices are equidistant to the nearest 
corner (i.e. where two coordinates are zero) of the grove. Groves are acyclic — every 
connected component of a grove is a tree. Lastly, each grove spans I[n). These 
connectivity properties are in fact what distinguish groves from arbitrary subgraphs 
of Q{n), and so give us a combinatorial definition of groves. 
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Figure 4. A grove g of order 5, superimposed on a picture of 1(5). 

Within a grove notice that there are two types of edges: long edges and short 
edges, depending on whether the long or short diagonal of a rhombus is used. For 
a vertex v = k) in C(n), we say that v is: 

up: if i + j + fc = — n + 1, 

down: if i + j + fc = — n — 1, 

fiat: \i i + i + k = — n, 

even: \i i+ j + k is even, 

odd: if z + j + A; is odd. 
Long edges connect flat vertices to flat vertices, and short edges connect up vertices 
to down vertices. Even vertices are only connected to even even vertices and odd 
vertices are only connected to odd vertices. It is shown in that every vertex in 
B{n) has degree 2 and only uses its short edges. As a result, there are only finitely 
many long edges, and these determine the grove. This observation leads to a more 
convenient way of looking at groves. 

1.2. Simplified groves. We begin by constructing a modified form of the cube 
recurrence. Let aij^ bkj, Ci^k be long edge variables where —n = i + j + k is fixed. 
The variable Uij is the label for the edge between vertices — l,fc + 1) and 
{i — k + 1), bkj is the label for the edge between {i — k + 1) and k), 
and Ci,k is the label for the edge between k) and — 1, fc + 1). We write a 
modified form of the cube recurrence as follows: 

fi,j^kfi—lj — l,k—l bi^J^Cijf-i^lJ^kfi.j — l^k—l ^ ei^jOjj^ f^j — ij^ f^—ij ^f^^—i 

+0,j,kbi.kfi,j.,k-lfi-l,j~l,k 

As we said, the long edges determine the grove, so rather than setting /ij,fc := Xij,k 
for {i,j,k) € we set fij.k ■= 1 for {i,j,k) S I{n). Then /o,o,o is simply 

a polynomial in the edge variables aij,bij,Cij, where the variables appear with 
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exponent +1 or 0, depending on whether the corresponding long edge is present or 
not. Each term describes a unique grove, and we stiU produce every grove. This 
form of the cube recurrence is caUed the edge variables version. 

Taking inspiration from the edge variables version of the cube recurrence, we 
can draw a simpler picture of our groves by ignoring all short edges and all of the 
vertices incident with them. In other words, specify a subset of the standard initial 
conditions of order n, called the simplified initial conditions: T' [n) = {(i,j,k) e 

\ i + j + k = —n,i, j,k < 0} CZ 1{n) . The simplified initial conditions are just all 
of the flat vertices. We now represent our groves as graphs on this vertex set — a 
triangular lattice shown in Figure |5l Also in Figure [S] we see the same grove as in 
Figure^ but with only the long edges included. In terms of edge variables, this 
grove is given by 

Go, 010400^211, 001402460, 0^*0,1 Co, oCo,lCl,oC2,0- 




Figure 5. On the left: X'(5) compared to X(5). On the right: a 
standard grove and its corresponding simplified grove. 

Another modification of the cube recurrence that we shall like to use is the 
edge-and-face variables version. In the original version of the cube recurrence, 
the variables Xij^k such that i + j + k = —n + 1 were vertex variables. In the 
simplified picture, we call them the face variables of order n, for reasons that 
will become clear. Rather than setting fij^k '■= 1 for all {i,j,k) in X(n), we give 
the face variables their formal weights. That is, we set fij^k '■— 1 for {i,j,k) e 
{{i,j,k) e I?\ - n - 1 < i+ j + k < n,i,j,k < 0} and /i,j,fc := Xi^^k for 
{i,j,k) e { (i, j, fc) e Z'^ I i + j + fc = —n + l,i,j,k < 0}. Generating /o,o,o using 
these initial conditions, we get a Laurent polynomial in the edge and face variables. 
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The vertices of the simphfied initial conditions can be seen as forming n(n+l)/2 
downward-pointing equilateral triangles, each with top-left vertex — 1, fc -I- 1), 
top-right vertex {i — l,j,k + 1), and bottom vertex {i,j,k). The face variables 
then correspond to each of these downward-pointing triangles. The triangle with 
{i,j,k) as its bottom vertex has face variable Xij_k+i- The exponent of the face 
variable is —1,0, or 1, corresponding to whether the downward-pointing triangle 
has, respectively, two, one, or zero edges present. There can't be three edges, since 
that would introduce a cycle and we would no longer have a forest. Although the 
face variables don't tell us anything new about a particular grove, they will be 
useful later in deriving probabilities of edges being present in random groves. 

1.3. Grove shuffling. We have given one definition for what groves are, and how 
they may be generated. The methods and notation introduced in the previous 
section will be very helpful for later proofs. However, there is another tool we 
will like to use; an algorithm called grove shujfling (or cube-popping as in W). 
Grove shuffling not only gives a purely combinatorial definition of groves, but also 
a method for generating groves of order n uniformly at random. Its inspiration 
comes from domino shujfling, due to Elkies, Kuperberg, Larsen, and Propp 
The use to which we put grove shuffling is directly motivated by James Propp and 
his paper JIB!. For proof that grove shuffling does indeed give rise to the same 
objects as the terms of the Laurent polynomials given by the cube recurrence, see 
Carroll and Speyer Here we will only include a description of the algorithm. 



Grove shuffling can be thought of as a local move on the downward-pointing 
triangles of a simplified grove according to whether a triangle has zero, one, or two 
edges present. See Figure |S1 Let a; be a generic downward-pointing triangle with 
possible edges a, b, c as shown, and let x' be an upward-pointing triangle, concentric 
with x, with possible edges a',b',c' as shown. There are three configurations of 





Figure 6. Grove shuffling. 
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X with two edges: ab, ac, be. Grove shuffling takes each of these triangles and 
replaces them with an upward-pointing triangle x' having none of its possible edges 
present. There are three configurations of x with exactly one edge: a, 6, c. Each 
of these is replaced by the upward-pointing triangle x' with only the parallel edge: 
a', 6', c', respectively present. Lastly, there is one configuration of x with none of its 
possible edges present. This triangle is replaced with the upward-pointing triangle 
x' containing any two of its three possible edges: a'b' ^a' c\b' c' , chosen randomly 
with probability 1/3. This last step is the only random part of the algorithm. After 
we have turned every downward-pointing triangle into an upward-pointing triangle, 
we add three new vertices to the corners of the grove so that we may shuffle again. 
For an example of grove shuffling, see Figured"^ 





Figure 7. A grove of order 4 shuffled into a grove of order 5. 

There is a unique grove of order 1. It has one downward-pointing triangle with 
zero edges. We now give a purely combinatorial description of simplified groves 
on standard initial conditions of order n: they are all the possible results of n — 1 
iterations of grove shuffling, beginning with the grove of order 1. From looking at 
the cube recurrence, it is not hard to show that there are 3L" /^J groves of order n. 
We can now make the following claim about grove shuffling. 

Theorem 1. Beginning with the unique grove of order one, any grove of order n 
will be generated after n — 1 iterations of grove shuffling with probability 1/3 L" /^J . 
In other words, grove shuffling can be used to generate groves uniformly at random. 

Proof. Clearly the statement holds for n—2. Suppose that the claim holds for some 
k > 1. We would like to know the probability of an arbitrary grove of order k -\- 1 
being generated. Fix such a grove and call it G{k -1-1). Only a certain subset of 
the groves of order k can be shuffled to become G{k -1-1). Call this set the shuffling 
pre-image of G{k + 1), denoted S-^{G{k + 1)). Let G{k) £ S-^{G{k + 1)). Let a 
be the number of downward-pointing triangles in G{k) with zero edges, let b be the 
number with exactly one edge, and c be the number of downward-pointing triangles 
with two edges. 

From the rules of grove shuffling, we see that the order of S^^{G{k + 1)) is 
3^. Each pre-image is obtained by making different choices of the the two edges 
appearing in each of the c downward-pointing triangles of G{k). So since we have 



■^To see grove shuffling in action, visit |http://ups. physics. wise, edu/ 'hal/SSL /groveshufflCTT] 
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supposed the probabihty of generating a particular grove of order k to be uniform, 
the probabihty is 

3L?T 

that after k shufELes we produce a grove in S~^{G{k + 1)). 

Let S{G{k)) = S{S~^{G{k + 1))) be the set of groves of order k + 1 that can be 
obtained by shuffling a grove in S-^{G{k + 1)). The order of S{G{k)) is 3". This 
is because in each of the pre-images there are a downward-pointing triangles with 
no edges present, and every such triangle can be shuffled to any of three upward- 
pointing triangles. Furthermore, the only edges where the groves of S~^{G{k + 1)) 
differ will be annihilated upon shuffling. So there is a 1/3° chance that one of the 
pre-images of G'(fc + 1) will actually shuffle into G{k + 1). Therefore the probability 
that k + 1 iterations of grove shuffling yields G{k -I- 1) is 

1 1 
SL^J ' 3^' 

Now we claim that a — c = [-^^J . If so, then the probability computed above is 
equal to 

1 



as desired. 

Let us make some basic observations from ^ or by easy induction. First, a + 
b + c = k(k + l)/2; the total number of downward-pointing triangles in any grove 
of order k. Secondly, b + 2c= [^J ; the total number of edges in any grove of order 
k. Then a — c = k{k + l)/2 — [^J ~ L^T^J i ^-i^d the theorem is proved. □ 

1.4. Frozen regions. We now describe the phenomenon that we analyze in Section 
121 First we observe that edges are indexed relative to the corners perpendicular to 
them, so in fact the edges a and a' in the description of grove shuffling have the 
same label: a = a' = aij. Horizontal edges are indexed relative to the bottom 
corner, and the diagonal edges are indexed relative to the top-right and top- left 
corners. In this way we can think of grove-shuffling as more akin to domino shuffling 
|lti|. Rather than replacing edges with parallel edges, we "slide" edges toward the 
corners along perpendicular lines. When a downward-pointing triangle has two 
edges, we remove both of those edges because they "annihilate" each other. When 
a downward-pointing triangle has no edges, we create two new ones randomly. 

With this viewpoint, we define an edge to be frozen if it cannot be annihilated 
under any further iterations of grove shuffling. Clearly the bottom corner edge, 
ao.o, is frozen when present. Then the edge aij is frozen exactly when the edges 
o>i' _j' are frozen, i < «' < 0, j < j' < 0. Diagonal edges behave similarly. In Figure 
ISlall the highlighted edges are frozen. 

We conclude this section by examining a picture of a large random grove gener- 
ated by grove shufHing. In Figure we see that outside of a certain region, all of 
the edges are parallel. Moreover, the boundary between the less uniform interior 
and the frozen regions in the corners seems to approximate a circle. Proving that 
this boundary approaches a circle in the limit is the main goal of this paper. 
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Figure 9. A grove on standard initial conditions of order 100. 



2. The arctic circle theorem 

For any n, we can scale the initial conditions so that they resemble an equilateral 
triangle with sides of length \/2, by mapping each vertex (z, j, k) to {i /n, j /n, k/n). 
The corner vertices (— n, 0, 0), (0, — n, 0), (0, 0, — n) are scaled to (—1, 0, 0), (0,-1, 0), 
and (0,0,-1). We will show that outside of the circle inscribed in this triangle, 
there is homogeneity of the edges in an appropriately scaled random grove of order 
n, with probability approaching 1 as n ^ oo. Specifically, we will examine the 
limiting probability of finding a particular type of edge in a given location outside 
of the inscribed circle. 

2.1. Edge probabilities. Let Pn{h j) = pihj, k), k = —n—i—j, be the probability 
that aij{n), the horizontal edge on triangle Xij^k+i, is present in a random grove 
of order n. Similarly define probabilities qn{k,i),rn{k,j) for the diagonal edges 
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bk,i{TL) and Ck.j{n) of the same triangle. Define En(i,j) — E{i,j,k + 1) = 1 — 
Pnihj) — (ln(k,i) — rn{k, j). The numbers En{i,j) are analogous to the creation 
rates discussed in [5], [2], and ^Hl- We will also refer to them as creation rates. As 
proven below, we can also realize the number E„{i,j) as the expected value of the 
exponent of the face variable Xi^^k+i- We prove the following formula for finding 
the edge probability Pn{hj) in terms of creation rates. ^ 



Theorem 2. The horizontal edge probabilities are given recursively by pn{i,j) 

2 

Pn-i{i,j) + |£'„_i(i,i). Thus, Pn{i,j) = -T^^EiiiJ). 



3 

1=1 



Proof. We wish to derive a relation between Pn{hj) and pn-i{i, j). In order to 
simplify notation, we let: 

P^Pn-i{i,j) q = qn-i{k,i) r = r„_i(fc,j) 

a = aij(n-l) b ^ bkain - I) c = Ckj{n-l) 

P^Pnihj) Q^qn{kJ) i? = r„(fc,j) 

A^aij{n) B = bk,i{n) C = Ckj{n) 

See Figure Uni 




Figure 10. Labels for downward- and upward-pointing triangles. 



Let pr{*), where * is a subset of {a, 6, c}, be the probability that a random grove 
contains that set of edges and not its compliment. Define Pr(*) similarly. Some 
observations that come directly from grove shuffling: 



priab) = pr{ac) — pr{bc) 

• priabc) — 

• pr{il}) + pr{a) + pr{b) + pr{c) + pr{ab) + pr{ac) + pr{bc) 
p — pr(a) -f pr{ab) + pr{ac) 
q — prib) + pr{ab) + pr{bc) 
r — pr{c) -f pr{ac) -f pr{bc) 
Pr{A) = pr{a) 
Pr\B) = pr{b) 
PriC) = prlc) 

Pr{AB) = Pr{AC) = Pr[BC) = l/3pr(0) 



^Notice the similarity between this statement and equation 1.5 of 



12 



T. K. PETERSEN AND D. SPEYER 



We will now deduce P = pn{i,j)- 

P = Pr{A) + Pr{AB) + Pr{AC) 

= pr(a) + 2/3pr(0) 

= pr(a) + 2/3(1 — pr{a) — pr{b) — pr{c) — pr{ab) — pr{ac) — pr{bc)) 

= pr{a) + 2/3(1 — p — q — r + pr{ab) + pr(ac) + pr(bc)) 

= pr{a) + 2/3{pr{ab) + pr{ac) + pr{bcj} + 2/3(1 — p — q — r) 

= pr{a) + pr{ab) + pr{ac) + 2/3(1 — p — q — r) 

= p + 2/5{l-p-q-r) 

Let X = be the face variable of the downward-pointing triangle in ques- 

tion. Notice that 

E{x) = Expected value of exponent on x 

= 1 • pr(0) + • {pr{a) + pr{b) + pr{c)) — 1 ■ {pr{ab) + pr{ac) + pr(bc)) 

= 1 — pr{a) — pr{b) — pr{c) — 2pr{ab) — 2pr{ac) — 2pr{bc) 

= 1 — p — q — r 

= E„_i{i,j). 

Therefore, P =p + 2/3E{x). In the coordinate system, we have 

2 "^^ 

Pn{i,j) = Pn-l{i,j) +2/3En-i{i,j) = ^i^^EiiiJ) 

1=1 

and the theorem is proved. □ 

2.2. A generating function. Wc now know that to compute the probability of a 
particular edge being present in a random grove, it will be enough to compute the 
creation rates Ei{i,j). In this section we derive a generating function for computing 
these numbers as well as the related generating function for the horizontal edge 
probabilities. 

Let F{x,y,z) = ^ E{—i,—j,—k)x'^y^z'' be the generating function for the 

i,j,k>0 

creation rates. First consider the uniformly weighted version of the cube recurrence: 

fi,j,k.fi-l,j-l,k-l = {,fi-lj,k,fij-l,k-l + fi,j-l,kfi-l,j,k-l + .fij,k-lfi-l,j-l,k) ■ 

We will return to the convention of setting fij^k = Xi,j,k for all (i, j, k) € I{n). Using 
this recurrence to calculate /o,o,o we will get each monomial weighted uniformly, so 
that if wc set all the variables equal to 1, /o,o.o = 1- If wc want the expectation of 
the exponent of the face variable x = Xi^^j^^ko, we need only calculate the derivative 
of /o,o,o with respect to this variable, then set all variables equal to one. In other 
words, 

E{io,jo,ko) = ^(/o,o,o) 

Furthermore, we can calculate the intermediate creation rates the same way. 

Lemma 1. Fix x = Xi„j„,k„ for io + jo + fco = — n + 1. Then for any {i',j',k') 
such that i' + j' + k' = —n' + 1 with n' < n, we have 

( fi' ,j' ,k'] = E{io - i',jo-j',ko - k'). 

ox \ / Xi,j,fe=l 



= 1 
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Proof. First, we re-center our initial conditions to make the situation clear. In- 
troduce the variables x by Xij^k — * Xi-i'j-j'^k-k'- In particular, Xi^jg^ko ~* 
Xio-i' ,jQ-j' M-k' ■ Then we have fi\j'M' /o,o,o, a polynomial generated from 
a set of standard initial conditions of order n — n' — 1. Differentiating fi'j^k' with 
respect to Xi^jg^ko and setting all the variables to one is equivalent to differentiating 
/o.o,o with respect to Xi^^ii ^ko-k' and setting all variables equal to one. The 
latter action clearly gives E{io — i' ,jo — j' ,ko — k'), and the claim is proved. □ 

With this in mind, let us differentiate the weighted cube recurrence with respect 

to X = ^io,jo)fco* 

fi,j,kfi-i,3-^,k-l + fi,j,kfi-ij-i^k-l = g (/i-1 j",fe/i,i-l,fe-l + /j-lj',fc./'jj-l,fe-l + 

fi,j-l,kfi~l,j,k-l + /jj-l,fe/j-lj,fe-l + 
fi,j,k-lfi-ij-i,k + fi,j,k-lfi-ij-i^k^- 

Now by setting Xij,k — 1 for all k), we get a linear recurrence for the expecta- 
tions in question (where r = io — i, s = jo — j,t = ko — k): 

E{r,s,t)+E{r + l,s + l,t+l) = ^(E{r + l,s,t) + E{r,s + l,t+l) + 

E{r,s + l,t) + E{r + l,s,t+l) + 
E{r,s,t+l) + E{r + l,s + l,t)y 

The recurrence holds for any r,s,t < 0. Let us also observe some recurrences near 
the boundary. 



E{r,0,0) 


= ^E{r + 1,0,0) 


E{0,s,0) 


= ^EiO,s + l,0) 


E{0,0,t) 


= ^E{0,0,t+1) 


E{r,s,0) 


= i (E{r + 1, s, 0) + E{r, s + 1,0) + E{r + l,s + 1, 0)] 


E{r,0,t) 


= ^(^E{r + l,0,t) + E{r,0,t + l) + E{r + l,0,t+l)'^ 


E{0,s,t) 


= i (^E{0, s+l,t)+ E{0, s,t+l)+ E{0, s + l,t+l)^ 



After computing £^(0, 0, 0) = 1, we can form the rational generating function in the 
variables x,y,z: 

F{x,y,z) = Yl E{-i,-j,-k)x'y^z>' 

i,j,k>0 

_ 1 

1 + xyz - l{x + y + z + xy + xz + yz) 
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Now using the fact that p{i,j, k) = p{i,j, fc + 1) + {2/3)E{i,j, k + 2), we can derive 
the formula for the generating function we want: 

G{x,y,z) = p{-i,-j,-k)x'y^z'' 

i,j,k>0 

2z^F(x,y,z) 
3(1 -z) 

2.3. Asymptotic analysis. With our generating function in hand, we can prove 
our main theorem. First let us embed a triangle in three-space by T := {{x, y, z) € 
E.^ \ x,y, z < 0,x + y + z = —1}. This is the triangle that we will scale to fit. 
A point {x, y,z) Q T is outside of the inscribed circle (what will show is the arctic 
circle) if and only if the angle between the vector {x, y, z) and vector (—1, —1, —1) 
is greater than cos~-^(y'2/3). 

Notice that for any point {x^y, z) outside of the inscribed circle, we have either 
x^y > —1/2, x,z > —1/2, or y,z > —1/2. We call any coordinate with a value 
strictly greater than —1/2 a small coordinate. We now state 

Theorem 3 (Weak Arctic Circle). Let (xo,yo, zq) be a point in T outside of the 
inscribed circle for which zq is a small coordinate. Let (i„, j„, kn), in-^jn+kn — —n, 
be a sequence of nonpositive integer triples such that 

lim -{injn,kn) = {xo,yo,zo) 

n—too n 

Then lim p{in,jn,kn) = 0. 

n — ^oo 

In other words, the theorem states that in the upper two regions of T outside of 
the arctic circle, the probability of finding a horizontal edge goes to zero as the order 
of a (scaled) random grove goes to infinity. By symmetry, there can be no diagonal 
edges in the lower region, and in order to satisfy the connectivity properties of 
groves, all the edges in the lower region must be horizontal. The following lemma 
is the heart of the proof. 

Lemma 2. Fix a point {xq, yo, zq) in T outside of the inscribed circle for which Zq 
is a small coordinate. Then there are real constants A, B, C such that 
p(-*,-J,-fc) = 0(e-(-4»+W)) 

for all i,j,k > and Axq + Byo + Czq < 0. 

Proving Lemma [21 consumes most of the rest of this section. Let us suppose the 
lemma is true and present the proof of the theorem. 

Proof of Theorem\^ By Lemma El p(i„, j„, /cn) = 0(e'^'"+-^-'"+'^''"), so we wih 
have that p{in,jn, kn) ^ if Ai„ + Bjn + Ckn —oo. 
Say Axq + Byo + Czq — d < and let 

a„ = A h B h C — . 

71 71 71 

Then for any e > there is some > such that for all ti > iV, a„ e B{d, e) = 
{x : \d — x\ < e}. So if we take n sufficiently large, na„ < 0, and \d + e|7i < |a„|7i. 
Since |d + e|n — > c», by comparison we have 

7ia„ = Ain + Bjn + Ckn -oo, 

and the theorem is proved. □ 
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To facilitate the proof of Lemma |21 we will use the following claims. 



Claim 1. Let f{x,y,z) be an analytic function. Let r,s,t be positive real numbers 
such that f{x, y,z) ^ for \x\ < r, \y\ < s, \z\ < t. If 



G{x,y,z) 



1 



f{x,y,z) 



i,j,k>0 

then aij,k — 0{r^^ s^H^^). 

Proof of Claim^ Define the loops 7 = {\z\ = t},7' = {|y| = s}, and 7" = {|x| = 
r}. Then we have 

1 d d d 



iljlkl dx^ dy^ dz^ 
1 



-G{x,y,z) 



< 



(27ri)3 

M 
(2^ 



7 "'7' "'7" 



(0,0,0) 
^^"^'^'^^ dxdydz 

„i+lyj + lzk+l "'^^ 



7 "'7' "'7' 



^y-i^^z'^ 
1 



2:i+lyi + l^fe + 



Y dxdydz 



= M 



(Since G{x, y, z) is bounded on the compact set 7 x 7' x 7".) 
1 



r* sH^ 

'i,j,k 



In other words, Oi ^ i. = 0{r 's H ^) and the claim is proved. 



□ 



Claim 2. Let {xQ,yQ, zq) be as in Lemma\^ Without loss of generality, say yo is 
the second small coordinate. Then there exists a vector {Aq, Bq,Cq) € R"^ such that 

• ^oa;o + Boyo + CqZq < 

• AoBq + AoCo + BoCo > 

• So, Co < 0. 

Proof of Claim\^ Let (j) be the angle between (cco, 2/0,^0) and (— 1,— 1,— 1). The 
assumption that {xo,yo,zo) lies outside the Arctic circle means > sm^^{y^l/3). 
Take (Ag, Bq^Cq) in the plane spanned by (— 1,— 1,— 1) and {xQ,yQ,zo) and ly- 
ing on the other side of (— 1,— 1,— 1) from {xo,yo,zo). Let 9 be the angle be- 
tween {Ao,Bo, Co) and (—1, —1, —1), so 9 + <j)is the angle between {Aq, Bq, Cq) and 
(a;o,yo,^o)- If = sm^^{^/^/3) -f e, e > 0, we can choose 9 = cos^-^(-\/l/3) — e/2 
so that 9 + 4) ~ cos 
From 9 < cos^^ 



^/l/3 + sm^^ yn/3 + e/2 > cos^i ^/l/3 + sm^ 
1/3, we deduce 

(^o + i3o + Co)^ ^ ((Ao,i3o, Co), (-1,-1,-1))^ 
3{Al + B^,+C',) 



A/3 = tt/2. 



||(Ao,i3o,Co)|P||(-l,-l,-l)|P 

^ cos^ 9 
1 

> 3 

so (Ao + Bo + Co)^ > A^ + B^ + C^ which is equivalent to ^o^o + ^oCo + BoCo > 0. 
From 9 + (j) > n/2, we deduce that AqXq + Boyo + CqZo < 0. 

Finally, we must show that we can choose Bo, Co < 0. Let a be the distance 
between (Ao,i3o,Co) and l/3(— 1, — 1, — 1). Let /3 > V3/6 be the distance between 
(xq, 2/0,2:0) and l/3(— 1, — 1, — 1). We have, for (Ao,i3o,Co) satisfying the first two 
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W I * 

i/ (-1/3,-1/3.-1/3) 





Figure 11. The relative positions of (xq, j/oj -zo) and (Aq, Bq) C'o)- 



conditions and lying in the plane x + y + z = —1, that q:/3 > 1/6. Now let 
u = {xq, yo, zq) + 1/3(1, 1, 1) so that \v\ = /3. Then choosing 

{Ao, Bo, Co) = -a^ - ^(1, 1, 1) = ~v- 1, 1) 

has all the desired properties. The first two are obvious by construction, and using 
the facts that —a < —1/(6/3), f3 > \/3/6 and, because zo is a small coordinate, 
—zo < 1/2, we can verify that 

--( 

-zo - 1/3 



Co 



< 



< 



{zQ + 1/3) - 1/3 
1/3 
1/3 



6/?2 
1/2-1/3 



6/3- 



1 



1/3 



36/32 
< 1/3-1/3 = 

If yo is the other small coordinate then an identical computation shows iJo < 0. □ 

Claim 3. Let {A, B, C) G with C < 0. Suppose that 

g{x, y, z) = 1 + xyz — (1/3) (a; + y + z + xy + xz-\' yz) 

is not zero for {x,y,z) G [0,6^^] x [0, e-^] x [0, e'-^]. Then g{x,y,z) is not zero for 
any {x, y, z) e with {x, y, z) e B(0,e^) x B(0,e^) x B(0,eC) = { {x, y, z) : \x\ < 
e^,|y|<e^,|^|<e^}. 
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Proof of Claim\^ Suppose for contradictfon that there is a complex zero {x,y,z) 
of g with Ij/I, Izj) S [0, e"*] x [0,e"^] x [0,6*-^], but no real zeros in the same 
region. As the zero locus of g is closed, we may assume that there is no complex 
zero {x',y',z') with < \x\, \y'\ < \y\ and \z'\ < \z\. So the power series of 

1 



G{x,y,z) 



(1 - z)g{x,y,z) 



converges on _B(0, |a;|) x B{0, \y\) x B{0, \z\) (we used C < to conclude that the 
(1 — z) term doesn't vanish) and blows up to cxi as we approach (a:, y, z). 

But the coefficients of G{x,y,z) are all positive as they are probabilities. So 
the series must also blow up as we approach (|a;|, \y\, \z\) and thus g{\x\, \y\, \z\) — 
contradicting our assumption that there are no zeroes in [0, e"*] x [0, e^] x [0, e*^] . □ 

Proof of Lemma. We now apply the claims to the edge probability generating func- 
tion: 

G{x,y,z) = ^ p{-i,-j,-k)x'y^z'' 

i,j,k>0 



(1 — z)(l + xyz — ^{x + y + z + xy + xz + yz)) 



2 

^2 



4(1 - z)g{x,y,z) 



2 > 

By ClaimHwe only need to show that we can choose real numbers A, B, C so that 

• Axq + Bya + Czq < 

• Both 1 — z and g{x, y, z) are not equal to zero for any (x, ?/, z) g { {x, y, z) e 
C3 : \x\ < e^,\y\ < e^,\z\ < e^}. 

We will now show that, for {Aq, Bf),Co) as in Claim [3 and t positive and suf- 
ficiently small, {A,B,C) = t{AQ, BqjCq) has the desired properties. We have 
Axq + ByQ + Czq < and C < because the analogous properties hold for 
{Aq, Bq, Co). All that remains is to show that g{x, y, z) does not vanish for (x, y, z) e 
B{0, e^) X -6(0, e^) x B{0, e*^). By Claim|21it is enough to show that g has no zeroes 
on [0,e^] X [0,6-^] x [O.e"^]. The identity 

(1 - x){l - yz) + y){l ~ xz) + (1 - z)(l ^ xy) 
g[x, y, z) ^ 

shows that g{x,y,z) 7^ on [0,1)'^. Writing x — e", y — e^, z — e'^ , we have 
g{x, y, z) = aP + a'^ + + 0{a + f3 + j)^. So, near (1, 1, 1), the zero locus of g 
looks like the cone a/3 -I- 07 -I- /?7 = 0. 

Let £ — {{x,y,z) : g{x,y,z) = 0} be the zero locus. We want to show that 
there is a number t > such that the point {e^, ,e'^) = (e*'*", e*^°, 6**^°) is inside 
of £. We can write 

ff(e*^°,e*^",e*^«) = t{AoBo + AoCo + BoCo) + t^O{Ao + Bq + C'o)^ 
= t{AoBo + AoCo + BoCo) + t^O{-l) 
= t{AoBo + AoCo + BoCo + t^O{~l)). 

For fixed {Aq, Bo,Cq), we can certainly choose a t > small enough to guarantee 
that AoBo + AoCo + BqCq + t^O{-l) > 0. Therefore we have (e^, e^, e^") ^ £ and 
moreover it is on the 'inside' of £ in that it is closer to the origin than the locus. 
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Now we'd like to say g will not vanish on [0,6"*] x [0, e^] x [0, e'-^]. Since it is 
nonzero on [0, 1)'^, we just need to check that g is not zero on [1, e^] x [0, e^] x [0, e*^], 
where we can take _B, C < 0. Let x = x{y, z) be a parameterization of the zero 
locus. Then we have 

_ l/3(y + z + j/z) - 1 

^^y'^^- 1/3(1 + y + 

It will be enough to show that for any (y, z) £ [0, e^] x [0, e*^], 

< x{e^,e^) < x{y,z). 

Let {y,z) be any pair such that < y < , < z < e'^ and < x{y,z). The 
pair (e^, e*^) is such an example. Then we will show that x{y, z ~ e) > x{y, z) for 
any < e < z, and similarly, x{y — e, z) > x{y, z) for any < e < We have 

l/3(y + z-e + y(z-e))-l 



x{y,z- e) 



y{z^e)-l/3{l + y + z-e) 
l/3{y + z + yz)-l-l/3e{y+l) 



yz - 1/3(1 + t/ + z)-e(y- 1/3) ' 

Because both y and z are less than one, the numerator and denominator of x{y, z) 
are both negative. It now suffices to check that e{y — 1/3) < l/'3e{y + 1). This 
amounts to saying y < 1, which is true by supposition. Hence, x{y, z — e) > x{y, z). 
By symmetry, x{y — e, z) > x{y, z) as well, and the lemma is proved. □ 

The asymptotics of multivariate generating functions is described in some gener- 
ality in the sequence of papers (21 1 Eli t)y Robin Pemantle and Mark Wilson. 
Unfortunately, their methods do not apply to our generating function. Refer to the 
complex 2-manifold given by 1 + xyz — l/3(x + y + z + xy + xz + yz) = as the 
singular variety: what we refer to above as the zero locus. Our analysis requires 
studying the behavior of this surface near the point (1, 1, 1). This point is a sin- 
gularity of the variety that, as we have seen, resembles a cone point locally. The 
papers of Pemantle and Wilson give asymptotics near smooth points of the singular 
variety and near multiple points of the singular variety, but not near singular points 
of the singular variety as in our case. We hope that an extension of their techniques 
can be used to prove a stronger version of Theorem In particular, we hope for 
a theorem that describes the asymptotic probabilities in any location (rather than 
just outside of the circle), similar to Theorem 1 of [2]. 

3. Domino tilings of Aztec diamonds 

We now draw parallels between the examination of the behavior of large groves 
on standard initial conditions, and the behavior of tilings of large Aztec diamonds. 
This approach yields no new results for Aztec diamonds, but presents an alterna- 
tive approach to their study. In this section we derive a generating function for the 
probabilities Pnih j) that position in a tiling of an Aztec diamond of order n is 
covered by a particular type of horizontal domino. Weak asymptotics for the func- 
tion we will derive are discussed as an example in )13j , and ongoing work of Henry 



^The introduction to 1131 refers to the case of a singular point of the singular variety as a 
tractable problem to be handled in future work. Personal communication from Pemantle reveals 
that work on the case where the singular point locally resembles a quadratic surface, e.g. as in 
our case of the cone, may be finished very soon. See section |21 and mention of current work of 
Cohn and Pemantle. 
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Cohn and Robin Pcmantlc seeks to give a full asymptotic expansion.^ The first 
derivation of the function is due to James Propp and Alexandra lonescu, though 
their (different) derivation has never been published. Some recursive formulas for 
Pn{i,j) are given in J^, and are the inspiration for our derivation of the edge 
probabilities for groves. 

3.1. Tilings of Aztec diamonds and the octahedron recurrence. We begin 
by describing precisely how tilings of Aztec diamonds are encoded in the terms 
of polynomials generated by the octahedron recurrence.^ First of all, rather than 
considering tilings of an Aztec diamond, we prefer to consider perfect matchings of 
the dual graph of this region. We call such a graph an Aztec diamond graph. The 
Aztec diamond graph of order n has as its vertices the set { (i ± 1/2, j ± 1 /2) | «, j S 
2j \i\ + |j| < n — 1 }. In other words, it is the set of centers of the unit squares 
that compose the Aztec diamond of order n. Each vertex is connected with an edge 
to its nearest horizontal and vertical neighbors. The faces of the Aztec diamond 
graph are the points (i, j) such that \i\ + \j\ < n. For example, see Figure [T^ 



(■1 



(P. 1.0) 









i-H— 



11-2 



Figure 12. Aztec diamond graphs of order n—1 and n — 2. 



Recall that the octahedron recurrence is: 

We initialize the octahedron recurrence by setting gij^n = Xi,3,n for n = 0,-1. 
Then 



tilings T 
of order n 



Thorough analysis of the generating function presented in this section requires analysis of the 
singular variety near a point that is locally a cone. 

^See Speyer's paper, 1171 . for more on encoding graphs with this recurrence. 
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where fn{T) is a Laurent monomial in the variables Xij,s where d = Oiii+j + n 
even, S = —1 ii i + j + n odd. Each monomial is of the form 

\i\ + \j\<n 

where deg(i, j) G {0, 1, 2} is the number of edges surrounding face (i, j) in the dual 
matching to the tiling T. For example, in Figure IT^ we see the tiling associated 
with the monomial 

3^-1,-1,03^1, -1,03^-1,1, 03^1,1,0 
2;o, -1,-12^0, 1,-12^0, 0,0 



Figure 13. Tiling of an Aztec diamond of order 2 and its dual matching. 

The assertion that Aztec diamond tilings can be so encoded is just a special case 
of Theorem 5.5 of Eric Kuo's paper on graphical condensation ^^l- Here we take 
the weight of a tiling T of order n to be m{T) as above. 

3.2. Domino shuffling. Domino shuffling is the name given to the local move that 
can be used to generate random tilings of Aztec diamonds. In a tiling of an Aztec 
diamond there are four types of dominoes, which may be characterized as follows. 
Make a checkerboard coloring of the Aztec diamond of order n by making the 
leftmost square in each row of the top half of the diamond be white. A horizontal 
domino is north-going if its leftmost square is white, south-going if its leftmost 
square is black. Similarly, a vertical domino is east-going if its topmost square is 
black, west-going if its topmost square is white. Domino shuffling is described in 
detail in |15) . but basically each domino slides in the direction indicated by its name. 
North-going dominoes take one step north, south-going dominoes head south, and 
so on. However, two dominoes may not pass through one another. If a north- and 
south- or east- and west-going pair of dominoes collide, they annihilate one another, 
and the resulting hole and any other holes opened by sliding dominoes is filled in 
with a pair of dominoes: two horizontals with probability 1/2, or two verticals with 
probability 1/2. Using this description we have. 

Theorem 4. The horizontal edge probabilities are given recursively by Pn{i,j) = 

^ n — 1 

Pn-i{i,j) + ^Er,-i{i,j). Thus, Pn{iJ) = - ^Ei{i,j). 

1=1 
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The proof of the theorem can be found in 16 . It foUows more or less dhectly 
from the definition of domino shuffling, where En{i,j) is the net creation rate (see 

0, nil). 

3.3. Another generating function. By differentiating the uniformly weighted 
version of the octahedron recurrence 

9iJ,n+l9iJ,n — l — ^ ^9i — ^,3,ngi+l,j,n + gi,j — l,ngi,j+l,n) , 

and because 
we obtain 

i(S„(^,J-l) + S„(^,J + l)). 
From this recurrence and Theorem 0] we get the generating function: 
G{x,y,z) = X! X! Pn{hj)x'y^z" 

">0|i| + |j|<« 

z/2 



(1 - yz){l + z2 - + x-i + y + y-^)) ' 

This is the form of the generating function used as an example in fOl. A weak 
arctic circle theorem like our Theorem|21follows directly from that example. Proba- 
bilities throughout the diamond could be extracted from this function in principle, 
and current work of Cohn and Pemantle seeks to carry out this more difficult anal- 
ysis. 

4. Biased groves, or, groves with a drift 




Figure 14. Groves with bias (a,/3,7). 
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Another variation on the model presented in this paper would be to study "groves 
with a drift," i.e., rather than having the random choice in grove shuffling be uni- 
form, we make a biased choice. * For example, say we choose the two diagonal edges 
with probability a, the horizontal and one of the diagonal edges with probability 
/3, and the horizontal and the other diagonal edge with probability 7 = 1 — a — /3. 
This bias is reflected in the cube recurrence by setting 

/ij,fc/i-l j-l,fe-l = Clfi-l,j,kfi,j-l,k-l + P fij -l,k fi~l,j ,k-l + 1 fij ,k-l fi-l,j -l,k- 

The generating function for biased creation rates is then 
Fa,i3.,',{x,y,z) = ^ Ea.i3,^{-i,-j,-k)x'y^ 

i,j,k>0 

__ 1 

1 + xyz — a{x + yz) — (3{y + xz) — 7(2 + xy) 
We can also derive the generating function for biased horizontal edge probabilities: 

Ga,0.n{x,y,z) = ^ pa,i3.,^{-i,-j,-k)x'^y^z'' 

i,j,k>0 

{l3 + j)z'^Fa^l3^^{x,y,z) 
1-z 

As seen in Figure^^ the arctic circle is just a special case of what one might call the 
"arctic ellipse," dependent on the free parameters a, (3. Using the same approach 
as in the unbiased case where a — (3 — 1/3, we can prove the following theorem 
without too much difficulty. 

Theorem 5. The boundary of the frozen region for (reseated) groves with a drift 
is given by the intersection of the plane x + y + z — —1 with the surface 

r'^ + s^ + <2 
rs + rt + st = ^ , 

where r — {f3 + "/)x, s — (a + 7)2/, and t = (a + l3)z. 

Aside from simply describing the shape of the frozen region in the biased situa- 
tion, it may also be useful to observe how the area of the temperate zone varies in 
a and /3 (with respect to the area of the entire grove) . Specifically, we can compute 
the ratio of the area of the temperate zone, Aq, to the area of the entire grove, A^, 
as a function of a and (3. Given that 7 = 1 — a — /3, 

f o^ ^(a + /3)(a + 7)(/3 + 7) 

((a + /3)(a + 7) + (« + + 7) + + l)iP + 1)) ' 

5. Speculation on statistics of groves 

As mentioned, we hope to apply the methods of Pemantle and Wilson to de- 
termine asymptotic probabilities throughout a random grove. Based on computer 
experiments and the similarity of groves and Aztec diamond tilings seen so far, 
we believe a formula for such probabilities exists. Another future aim is to apply 
the methods of growth models and statistical mechanics to groves, in the style of 
Johansson [7], or more recently Kenyon, Okounkov, and Sheffield |H1, 0) [Tn| . 



^See |5] for analysis of the same variation in Aztec diamond tilings. 
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5.1. Randomly growing Young diagrams. Perhaps something can be proved 
about the variance of the boundary circle by interpreting groves in a more familiar 
setting. One clever way for determining the boundary of the frozen region for Aztec 
diamond tilings is to look at a frozen corner as a randomly growing Young diagram. 
See [5], [H], [7j for this interpretation. A nearly identical projection of the frozen 
region of a grove yields some sort of randomly growing Young diagram, but it seems 
to follow more intricate rules of growth than those of Aztec diamond tilings. 



If we project the grove onto the plane ]R'^/(0,0, 1) we see a triangular array of 
boxes with two types of diagonal edges. Let us put the corner box at the upper-left 
and index it as (0,0). Then we index the rows by < i < n — 1, the columns by 
< j < n — 1, so that each box has index («, j) with i + j < n ~ 1. A diagonal 
edge from the top left corner to the bottom right corner of a box corresponds to 
a short edge in the grove. Edges from the top right to bottom left correspond to 
long edges in the grove. The box {i,j) is frozen if it contains a long edge and all 
the boxes contain long edges, < i' < i, < j' < j. Clearly the collection 

of all frozen boxes is a Young diagram. We would like to be able to describe how 
this Young diagram grows under grove shuffling. 

In a randomly growing Young diagram, we call box a growth position if 

boxes (i — l,j) and {i,j — 1) are both frozen and {i,j) is not frozen (we use the 
convention that boxes of the form (i, —1), (— 1, j) are always frozen). In the case 
of Aztec diamond tilings, the growth of a frozen corner under domino shuffling 
corresponds exactly to adding a new box at each growth position independently 
with probability 1/2. With groves this is not the case, though it may not be clear 
from Figure El In fact, two groves may have the same Young diagram projec- 
tion, yet grow very differently upon shuffling. It can happen that, at a particular 
growth position, one grove will not permit the addition of a new box under any 
circumstances, whereas the other will add a box with probability 2/3. 

Another difference seems to suggest the need for a new definition of growth 
position. It is possible that a grove can project to a Young diagram with growth 
position {i,j), but after just one iteration of grove shuffling the projection adds a box 
not only to position (i, j), but also position {i+l,j). In fact for any k there is a grove 
of some (perhaps large) order n with growth position (i, j), so that with positive 






Figure 15. Projecting a frozen corner to a Young diagram. 
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Figure 16. On the left, the Young diagram corresponding to a 
random grove of order 100. On the right, the Young diagram 
corresponding to a random tihng of an Aztec diamond of order 
100. 

probabihty its projection to a Young diagram adds boxes (i, j), (j+1, j), . . . , (j+fc, j) 
upon one iteration of grove shuffling. Perhaps such situations are outliers, but there 
is still much work to do in this direction. 

5.2. The nexus. Much of the statistical study of groves is motivated by analogy 
with statistics for domino tilings of Aztec diamonds. But there is at least one 
interesting feature of groves that seems to have no analogy in the realm of Aztec 
diamonds. We conclude the paper with some observations about a unique vertex 
that is present in every grove. We call this special vertex the nexus. Loosely 
speaking, the nexus is the vertex at the "middle" of the unique tree connected 
to all three sides of the initial conditions. If the nexus and its incident edges are 
removed from the grove, then all three sides become disconnected from one another. 
In Figure El we see the nexus highlighted. 




Figure 17. The nexus of a grove. 

Definition 2. Let g be a grove of order 2n or 2n+l, n = 1,2, .. .. The nexus of g is 
the (even) vertex v which is connected to each of the three midpoints.' (— n, — n,0), 
(— n,0, — n), and (0,—n,—n), and for which each midpoint lies on a distinct branch 
of the tree rooted at v. 

We would like to understand how the nexus moves during grove shuffling. The 
nexus takes some kind of random walk in the initial conditions, but it is not a simple 



AN ARCTIC CIRCLE THEOREM FOR GROVES 



25 



random walk. In some sense the nexus is a "stuttering" random walker. In some 
situations the nexus takes a step in one of three directions with equal probability, 
in others it does not move at all, and in still others it moves deterministically. 
Specifically, if the nexus is at an "up" vertex (see section []~T|l then after one iteration 
of grove shuffling it will take a step in one of three directions, each with equal 
probability. If the nexus is at a "down" vertex, then after one iteration of grove 
shuffling the nexus will not move (and so becomes a "flat" vertex in the new grove). 
Strangely, if the nexus is at a flat vertex then depending on which edges lead from 
the nexus to the midpoints, the nexus will either move deterministically one step 
or it will not move. Ultimately we want to say something about the distance of the 
nexus from the center of the grove as the size of the grove gets very large. How far 
from home does the nexus roam? 
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